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ABSTRACT 


Several  manpower  planning  models  are  presented  that  exploit 
the  longitudinal  stability  of  manpower  cohorts.  The  manpower 
planning  process  is  described  along  with  the  problem  of  identi¬ 
fying  and  obtaining  various  types  of  longitudinal  data.  An 
infinite  horizon  linear  program  for  calculating  minimum  cost 
manpower  input  plans  is  found  to  have  a  straightforward  solution 
in  a  great  many  cases  and  to  yield  an  easily  implemented  approxi¬ 
mation  technique  in  other  cases. 


1.  INTRODUCTION 


This  paper  formulates  manpower  planning  models  for  a  system  consisting  of 
many  skill  categories;  the  particular  application  is  to  the  enlisted  force  in 
the  U.S.  Navy.  Interactive  computer  models  of  these  theoretical  formulations 
have  been  developed  and  implemented  to  aid  decision-makers  who  wish  to  test 
the  effects  of  alternative  manpower  policies  on  staffing  requirements  and  fu¬ 
ture  manpower  budgets.  These  interactive  manpower  planning  models  have  a 
variety  of  uses,  to: 

(1)  predict  the  manpower  requirements  that  will  be  fulfilled  by  the 
current  stock  of  manpower. 

(2)  calculate  unfulfilled  requirements  and  the  new  Inputs  necessary  to 
meet  them. 

(3)  identify  bottlenecks  in  the  manpower  planning  process. 

(A)  assist  in  preparation  of  future  manpower  budgets. 

(5)  simulate  the  effects  of  manpower  policy  changes  on  future  manpower 
needs. 

(6)  relate  alternate  personnel  retention  and  performance  assumptions  to 
the  need  for  future  inputs. 

(7)  calculate  minimum  cost  input  schedules  when  lower  bounds  on  require¬ 
ments  are  given. 

Details  of  the  interactive  positions  of  the  models  are  not  included  here  since 
they  are  to  a  large  extent  dependent  on  the  particular  hardware/software  com¬ 
bination  used.  Rather  we  emphasize  the  underlying  mathematical  structure. 

Individuals  in  any  skill  category  can  be  identified  by  several  character¬ 
istics:  examples  are;  rank,  salary,  number  of  years  of  experience  in  the  skill 
category,  length  of  service  in  the  Navy,  and  personal  attributes  such  as  age 
and  measures  of  performance.  The  models  present  -d  in  this  paper  are  designed 
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to  assist  in  preparing  manpower  budgets  and  meeting  aggregate  strength  require¬ 
ments.  For  these  purposes  we  have  chosen  to  identify  individuals  in  the  enlisted 
force  according  to  skill  category  and  length  of  service  (LOS)  in  the  Navy. 

Section  2  of  the  paper  describes  the  underlying  manpower  flow  models. 

The  models  are  based  on  the  assumption  of  longitudinal  stability  in  the  service 
lifetimes  of  different  manpower  cohorts.  We  show  that  the  accession  schedule 
that  exactly  meets  requirements  is  found  by  solving  a  set  of  lower  triangular 
system  of  linear  equations.  Section  3  relates  several  methods  of  describing 
the  longitudinal  behavior  of  manpower  cohorts  and  shows  how  the  flow  model's 
parameters  can  be  estimated  from  existing  Navy  data.  In  Section  4  we  present 
an  in finite  horizon  linear  program  for  the  calculation  of  future  manpower  in¬ 
puts  and  strength  levels  with  the  objective  of  minimizing  discounted  costs. 

We  derive  readily  verifiable  conditions  on  the  inputs  to  the  infinite  horizon 
problem  that  guarantee  that  the  equality  solution,  described  in  Section  2,  will 
be  optimal.  In  cases  where  the  equality  solution  is  not  optimal  we  describe 
a  finite  linear  program  that  will  approximate  solutions  of  the  infinite  horizon 
program. 

The  models  presented  in  this  paper  examine  the  relationships  between  three 
factors:  (i)  the  current  manpower  situation  as  described  by  the  LOS  distribu¬ 
tions  of  different  skill  ratings,  (ii)  the  survivor  fractions  that  will  deter¬ 
mine  the  longitudinal  behavior  of  manpower  cohorts,  and  (iii)  the  manpower 
requirements  for  future  times.  The  size  of  our  models  and  the  type  of  calcu¬ 
lations  performed  allow  policy  makers  to  quickly  analyze  the  impact  of  various 
assumptions  and  policies. 


2.  THE  UNDERLYING  COHORT  FLOW  MODELS 

We  consider  an  organization  which  is  divided  into  many  skill  categories 
„i  =>.re  movement  between  categories  typically  involves  retraining  of  the  indi¬ 
vidual.  In  this  paper  we  consider  each  skill  category  separately;  a  subsequent 
paper,  will  di  ?ciuss  transfers  and  other  interactions  between  skill  cate¬ 
gories. 

We  idealize  the  evolution  of  the  skill  category  by  analyzing  its  changes 
at  discrete  points  in  time  (i  ■  ...  -2, -1,0, 1,2,  ...)  .  We  say  that  period 
i  is  the  interval  between  times  (i  -  1)  and  i  ;  it  is  a  future  period  if 
i  >  1  ,  a  past  period  if  i  <  0  ,  and  the  current  period  if  i  **  1  .  In  period 
i  a  number  of  people  enter  the  skill  category;  that  group  is  called 

cohort  1  and  x^  is  the  number  of  accessions  in  cohort  1  .  Let  be 

the  fraction  of  the  cohort  entering  in  period  i  which  is  still  present  and 
available  to  meet  requirements  at  time  i  +  j  (j  >_  0)  .  Let  z^  be  the  re¬ 
quirement  in  the  skill  category  at  time  k  and  let  (m  +  l)  be  the  maximum 
number  of  periods  an  individual  is  allowed  to  stay  in  the  system.  Thus 
=0  if  j  >  m  .  For  some  future  time  k  we  have 


(1) 


'k  *  Vk,0  +  ViVi.i 


+  . . .  + 


Xk-mak-m,m 


Equation  (1)  simply  says  that  the  requirement  at  time  k  is  made  uj  of  frac¬ 
tions  of  cohorts  which  survive  from  earlier  periods.  Thus  it  is  natural  to 
call  the  n£j,s  t^e  survivor  fractions  for  the  cohort  which  enters  at  time  i  . 
At  time  0  ,  the  history  of  past  accessions  is  given  by  the  vector 

(x  ,x,  ,  ...»  x  _,x  ,  ,xA)  .  The  current  inventory  of  people  is  given  by 

— m  x-ra  —4  —l  u 

x.a  -  +  x,a,,  +  ...+x  a  ,  and  contains  the  remaining  fractions  oi 
(J  U,U  “i  "l) X  *"in  in y in 

past  inputs.  This  quantity  is  called  the  current  legacy ,  say  y^  .  In  future 
period  k  the  legacy  y^  from  past  inputs  up  to  and  including  period  0  will 
be 


Jfrri. . 
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y,  ■  xnan  ,+x,a,,.,  +  ...+xf  a  ,  ,  if  k  <  m  , 

Jk  0  0,k  -1  -l,k+l  k-m,  -m,k+m  - 


if  k  >  m  . 


Suppose  we  have  a  planning  horizon  of  T  periods  with  requirements 
z^z^,  •••  2^  •  From  equations  (1)  and  (2)  we  see  that  future  cohort  sizes 
must  satisfy 


°l,0xl 


=  Z1  '  yl 


°1,1X1  +  °2,0X2 


Z2  '  y2 


“l,T-lXl  +  °2,T-2X2  +  •**  +  “T,0XT  "  ZT  “  yT  * 


Here  we  have  assumed 


Al:  requirements  are  met  exactly. 


Under  Al  it  is  quite  possible  that  for  a  given  set  of  z^'s  »  y^'8  an<* 
Cluj's  some  x^  could  be  negative.  Such  a  result  would  say  that  in  order 
to  exactly  meet  requirements  in  all  periods  1,2,  ...,  T  it  will  be  necessary 
to  remove  people  from  the  skill  category  in  period  k  .  In  practice,  this 
might  be  accomplished  through  retraining. 

This  section  concentrates  on  the  equality  solution  for  several  reasons. 
First,  it  is  misleading  to  state  the  problem  as  if  the  accessions  (the  x^) 
are  the  only  variables  which  the  decision-maker  can  influence.  The  legacies 
(y^)  ,  requirements  (z^)  and  to  some  extent  the  survivor  fractions  (a  ) 


Recall  that  a. 


0  if  j  >  m  . 
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can  all  be  changed  or  explicitly  influenced  by  manpower  policies.  Second, 
plans  that  are  eventually  recommended  will  probably  conform  to  the  equality 
constraints  since  budget  restrictions  do  not  generally  allow  for  slack  in  the 
system.  Third,  we  intend  to  use  the  models  in  this  section  to  test  the  ef¬ 
fects  of  alternate  policies  on  several  objectives:  (1)  the  departure  of 
realistic  requirements  from  ideal  requirements,  (2)  the  smoothing  of  stocks 
and  accessions,  (3)  the  impact  of  policy  changes  on  personnel  and  finally 
(4)  the  retraining  costs  associated  with  switching  personnel  from  one  spe¬ 
cialty  to  another.  In  Section  4  we  shall  drop  A1  and  treat  the  net  require¬ 
ments  as  lower  bounds  and  look  for  cost  minimizing  accession  schedules. 

In  the  remainder  of  this  paper  we  make  an  important  second  assumption 


A2:  the  survivor  fractions  a  are  stationary  from  period  to  period. 

i>j 

That  is,  a  ~  a  t  independent  of  i  and  independent  of  x.  . 

J  J  i 

Under  assumption  A2  equation  (3)  to  simplify  to 


U3X1 


(4) 


Vl  +  “ox2 


Vlxl  +  “T-2X2  +  •••  +  Vt  ■  ZT  '  yT 


In  these  equations  the  legacies  are  given  by 


yi  =  Vo  +  V-i  + 


+  a  x. 
m  1-m 


y2  =  Vo  +  V-l  + 


+  a  x 


m  2-m 


yT  = 


aTX0  + 


a  x„ 
m  T-ra 


(5) 
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Equation  (4)  ran  be  used  in  a  number  of  ways.  We  have  mentioned  already 
that  given  the  requirements,  legacies  and  survivor  fractions,  (4)  can  be  used 
to  calculate  new  cohort  input  requirements  for  each  period  of  the  planning 
horizon  T  .  Alternatively,  given  planned  inputs  over  the  next  T  periods 
the  z^s  can  be  considered  as  the  result  of  these  inputs.  Also,  given  re¬ 
quirements  and  planned  inputs,  the  legacies  which  satisfy  equation  (4)  can 
be  determined.  These  possibilities  will  be  discussed  further  in  the  section 
on  interactive  programming. 

To  this  point  we  have  motivated  equations  (4)  and  (5)  by  considering  z . 
and  y^  as  stocks  of  people  and  x^  as  a  flow  of  people  per  year.  It  is  not 
necessary  to  define  o  ,  x  ,  y  ,  and  z  in  that  manner.  We  could  also  speak 
of  stocks  and  flows  of  money  or  effective  personnel  or  use  another  unit  of 
time.  For  example  if  c^  is  the  cost  associated  with  an  individual  in  the 
j  year  of  service,  the  total  cost  in  year  i  will  be 


c ftoinx.  +  c,a,x.  -  ...  +  c  a  x. 

0  0  i  11  i-1  m  m  i-m 


The  cost  legacy,  that  will  be  incurred  in  future  period  i  due  to  current, 
period  0  ,  manpower  legacies  is 


(7)  c.a.xn  +  c.J,a.,,x.I  ...  +  c  a  x. 

i  i  0  i+1  i+1  1  mm  i-m 

The  use  of  equations  (4)  and  (5)  depends  on  the  availability  of  data  to 
estimate  the  various  parameters.  The  next  section  examines  several  possible 
data  configurations  and  presents  some  actual  examples  of  U.S.  Navy  enlisted 


force  data 
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3.  DATA 

To  employ  the  cohort  model  described  in  Section  2,  it  is  necessary  to 
obtain  information  about  past  accessions  x^  and  some  knowledge  about  past 
survivor  fractions  .  Typically,  data  which  is  available  in  most  organi¬ 
zations  is  of  the  cross-sectional  type.  Thus  we  may  know  the  number  currently 
in  a  given  skill  category  and  the  length  of  time  each  person  has  been  in  the 
category.  In  the  case  of  the  U.S.  Navy  enlisted  force  this  type  of  cross- 
sectional  information  is  available  for  a  number  of  past  planning  periods. 

Data  on  the  past  inputs  and  the  survivor  fractions  are  not  available,  end 
so  ve  must  either  eliminate  these  from  our  equations  or  find  methods  of  de¬ 
termining  them  from  available  data. 

We  shall  demonstrate  below  how  estimates  of  the  survivor  fractions  and 
past  skill  category  accessions  can  be  obtained  from  several  different  config¬ 
urations,  and  present  an  example  using  U.S.  Navy  data.  The  reader  should, 
however,  not  focus  too  strongly  on  the  problem  of  estimating  survivor  frac¬ 
tions  from  historical  data.  Survivor  fractions  realized  by  past  cohorts  are 
the  result  of  former  manpower  policies,  and  the  cohort's  behavioral  charac¬ 
teristics.  In  the  context  of  the  U.S.  Navy  these  data  are  influenced  by  early 
out  policies,  retraining  and  transfer  decisions,  changes  in  the  length  of 
tours,  the  mixture  of  regular  and  reserve  forces,  and  the  fraction  of  draft 
motivated  accessions.  Simple  projections  using  historical  estimates  of  sur¬ 
vivor  fractions  will,  at  best,  show  the  results  of  continuing  past  manpower 
policies.  The  specification  of  future  survivor  fractions  is,  to  a  great  ex¬ 
tent,  the  specification  of  future  manpower  policy. 

Despite  this  warning  on  the  use  of  historical  estimates,  we  shall  now 
show  how  these  estimates  can  be  made  with  several  different  data  configura¬ 
tions.  In  the  first  case  we  assume  that  past  accession  levels  x^  for  i  0 
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are  known  as  well  as  part  stocks 
with  assumption  A2  can  be  written 


Equation  (1)  of  Section  1,  combined 


(1)  2J  m  VO  +  Vl“l  +  XJ-2“2  •••  +  • 

If  2^  is  known  for  j  ■  0,-1,  -m  and  for  j  =  0,-1,  . ..,  -2m  , 

then  system  (1)  represents  m  +  1  equations  in  the  m  +  1  unknowns 
a0'al*  am  *  However»  is  highly  unlikely  that  data  on  skill  category 

inputs  will  be  maintained  at  the  precise  degree  of  aggregation  called  for  by 
the  model  over  the  past  2m  periods.  We  might  find  the  aggregate  input  into 
several  skill  levels,  or  a  change  in  manpower  accounting  methods  in  the  past 
2m  periods.  Thus,  in  general,  it  is  necessary  to  estimate  the  survivor  frac¬ 
tions  without  knowledge  of  past  accessions. 

We  are  now  left  with  the  task  of  using  the  models  described  in  Section  2 
when  past  accessions  are  not  known  and  cross-sectional  data  is  available.  We 
define  a  to  be  the  number  of  people  in  year  i  with  j  years  of  service. 
From  the  cohort  assumption  A2, 


(2) 


ij 


An  individual  has  j  years  of  service  if  his  length  of  service  (LOS)  lies  in 
the  interval  [j,j  +  1)  .  It  does  not  follow  that  an  individual  with  LOS 
equal  to  j  has  been  in  a  skill  category  for  j  years.  It  is  possible  to 
be  unrated  for  several  years,  and  it  is  possible  (o  be  transferred  between 
skill  categories.  In  many  of  the  critical  highly  trained  ratings,  however, 
the  LOS  and  time  in  skill  category  coincide. 

Note  that  formula  (2)  allows  one  to  calculate  the  legacy  without  knowledge 
of  the  past  accession  levels.  Given  th?  current  longitudinal  data  ngj  and 
the  survivor  fractions  that  will  obtain  in  future  periods,  formula  (2) 

can  be  used  to  calculate  the  past  accession  levels  that  would  be  consistent 
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with  both  current  LOS  data  and  the  assumption  that  the  survivor  fractions 
applied  in  past  periods.  These  accession  levels  can  then  be  used  in  equation 
(5)  of  Section  2  to  calculate  the  legacy. 

It  remains  to  find  some  method  of  estimating  the  survivor  fractions  and 
of  relating  them  to  other  measures  of  longitudinal  stability  namely:  the  con¬ 
tinuation  rate  0^  ,  and  the  age  fraction  y^  .  These  measures  will  be  de¬ 
fined,  and  then  related  to  the  survivor  fractions  ,  and  cross-sectional 
data  n^j  . 

The  continuation  rate  0.  is  defined  as 


n 


ill 


H,j-1 


It  is  most  convenient  to  interpret  0^  as  the  fraction  of  those  with  j  -  1 
years  of  service  in  year  i  -  1  ,  that  remain  (continue)  in  year  i  .  It  is 
possible,  although  this  effect  is  not  explicitly  considered  in  our  model,  to 
enter  a  skill  category  with  more  th'n  one  year  of  service.  Therefore  histori¬ 
cal  estimates  of  0^  from  reliable  data  can  yield  numbers  that  are  greater 
than  one.  The  notation  has  anticipated  the  result  that  0j  is  independent 
of  i  due  to  assumption  A2 .  From  (2) ,  we  have 


(4) 


Xi-ial 


i-l-(j-l)aj-l  aj-l 


j  =  1,2,  ...,  m 


The  value  of  0q  *  1  is  assumed.  Thus  we  can  calculate  0  from  o  using 
(4)  and  the  reverse  calculation  is  simply: 


(5) 


j  =  1,2, 


m  . 


The  age  fraction  y^  is  the  fraction  of  people  in  year  i  that  have  j 
years  of  service.  It  follows  that 


lO 


'iJ 


where  z.  a  £  n,,  .  When  assumption  A3  holds,  then 

j*0  13 


,  Xi-J  aj 
rij  z*  ‘ 


To  be  a  useful  measure  of  longitudinal  stability,  the  y^  should  be  sta¬ 
tionary,  i.e.  independent  of  i  .  However,  stationary  age  fractions  are  in 
general  incompatible  with  assumption  A2.  In  particular,  if  A2  holds  and  y^ 
is  independent  of  i  ,  then  one  can  easily  determine  that  the  system  must  be 
growing  at  a  constant  rate  (0  -  1)  with  ,  z ^  «*  G^Zq  ,  and 


13  0^ 


We  have  seen  that  the  survivor  fractions  a  ,  and  continuation  rates  3  , 
are  useful  measures  of  longitudinal  cohort  stability,  while  the  age  fraction 
does  not  isolate  the  effects  of  varying  past  levels  of  accessions  from 
the  behavior  patterns  of  the  cohorts.  We  conclude  this  discussion  with  a 
formula  that  shows  how  the  quantities  a  ,  3  ,  y  are  calculated  from  the 
cross-sectional  data.  Since  the  a  and  3  are  essentially  flow  rates,  at 
least  two  years  of  cross-sectional  data  are  necessary  to  determine  their  values. 


j  Vl,j-1  *i-lYi-l,j-l 


ni,1  ni2*nil _ o 

Vl.j-l  **•  ni-l,l'ni-l,0  0  ’ 


In  what  follows  we  shall  demonstrate  how  these  concepts  can  be  applied  to 
some  actual  data. 


11 


The  flow  models  described  in  Section  1  were  designed  to  assist  in  man¬ 
power  planning  for  the  U.S.  Navy's  enlisted  ranks.  That  force  comprises 
approximately  550,000  people  in  nine  ranks  and  £.0  skill  categories  which 
are  commonly  called  ratings. 

Tables  1  and  2  contain  the  length  of  service  distributions  of  five  ratings. 
If  we  take  June  30,  1972  to  be  time  0  ,  then  Table  1  gives  the  values  of 
n_^  j  and  Table  2  the  values  of  n^  ^  .  The  particular  data  presented  in 
Tables  1  and  2  was  abstracted  from  Rating  Career  Factor  reports  prepared  by 
the  Personnel  System  Research  Department  of  the  Naval  Personnel  Training  and 
Research  Laboratory,  San  Diego. 

The  Rating  Career  Factor  reports  contain  the  tail  of  the  age  fraction  as 
well  as  the  total  inventory  z ^  in  each  rank.  Let  T ^  be  the  fraction  in 
period  i  ,  that  have  completed  j  or  more  periods  of  service.  Clearly 

ri,j  "  Yi,j  +  YU+1  •**  +  Yi,m 

and 


y  =  r  _  r 
Yij  ‘i,j  i,J+l 


Thus  to  find  n^  we  must  calculate 


nij  ■  2i(ri,j  -  ri,j+i> 


For  the  Boatswain's  Mate  rating  we  have 


n0,5  =  z-l^r-l,5  “  r-l,6^ 

=  (10973) (5.97  -  56.2)  =  384  . 


The  original  source  of  this  data  is  the  enlisted  personnel  master  tape, 
that  is  maintained  by  the  Management  Information  Division  of  the  Bureau  of 
Naval  Personnel  (Pers  N) .  The  data  kept  by  Pers  N  is  used  on  a  daily  basis 
for  current  management  of  the  Navy;  it  is  a  necessity  of  the  cross-sectional 
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TABLE  1 

LENGTH  OF  SERVICE  DISTRIBUTIONS  AS  OF  6/30/71 


Length  of  Service 

BM 

QM 

Ratings 

ET 

CS 

I 

BT 

I 

0 

225 

1186 

5044 

1435 

1488 

1 

1547 

1078 

3848 

1372 

1926 

2 

1513 

911 

3544 

1173 

2156 

3 

382 

82 

2305 

357 

198 

4 

438 

141 

1762 

503 

285 

5 

449 

87 

‘  762 

493 

340 

6 

438 

125 

653 

399 

318 

7 

449 

141 

544 

325 

285 

8 

561 

152 

479 

378 

285 

340 

9 

516 

130 

436 

315 

10 

449 

168 

414 

284 

274 

11 

404 

125 

371 

231 

329 

12 

348 

130 

349 

231 

351 

13 

326 

130 

262 

273 

351 

14 

471 

168 

197 

451 

548 

15 

516 

168 

218 

545 

559 

209 

16 

236 

103 

.131 

294 

17 

427 

93 

110 

503 

307 

18 

483 

103 

88 

420 

220 

19 

124 

28 

44 

116 

1 

56 

20 

57 

17 

23 

32 

23 

21 

102 

17 

23 

43 

23 

22 

158 

28 

44 

64 

23 

23 

102 

12 

23 

32 

1 

24 

124 

17 

23 

64 

23 

25 

79 

17 

23 

32 

23 

26 

91 

17 

23 

32 

12 

27 

102 

17 

23 

53 

1 

28 

79 

17 

1 

32 

12 

29 

35 

6 

1 

11 

1 
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TABLE  2 

LENGTH  OF  SERVICE  DISTRIBUTIONS  AS  OF  6/30/72 


Length  of  Service  BM 

0  166 

1  769 

2  2832 

3  319 

4  341 

5  385 

6  418 

7  418 

8  418 

9  484 

10  495 

11  429 

12  395 

13  330 

14  319 

15  462 

16  495 

17  231 

18  297 

19  253 

20  78 

21  45 

22  67 

23  122 

24  78 

25  100 

26  67 

27  78 

28  89 

29  23 


Ratings 


QM 

ET 

CS 

BT 

1317 

3570 

1649 

2501 

596 

3578 

830 

961 

1050 

3578 

1489 

2026 

107 

2996 

331 

286 

62 

2091 

321 

170 

132 

1121 

470 

254 

87 

647 

480 

317 

117 

518 

361 

275 

122 

454 

301 

254 

137 

432 

351 

254 

127 

410 

311 

317 

157 

389 

291 

254 

117 

367 

231 

317 

132 

324 

231 

339 

127 

238 

261 

349 

167 

195 

441 

539 

167 

216 

530 

560 

97 

130 

291 

201 

62 

87 

311 

170 

56 

23 

191 

Ob 

16 

23 

71 

33 

11 

23 

21 

12 

16 

23 

31 

22 

21 

44 

41 

12 

11 

23 

31 

12 

6 

23 

41 

12 

16 

23 

21 

22 

11 

1 

41 

1 

21 

23 

41 

12 

6 

1 

11 

1 

type.  That  is,  on  any  given  day,  their  data  file  gives  a  current  listing  of 


Navy  Personnel.  We  obtained  values  of  n  .  .  from  report  E1S-9H,  Pers  .M2, 
Bureau  of  Naval  Personnel,  dated  9/30/71.  These  values  conformed  closely  to 
the  n  ,  .  calculated  from  the  age  tails.  From  Tables  1  and  2,  we  see  that 
on  June  30,  1971  there  were  1513  people  in  the  category  "Boatswain's  Mate"  in 
the  Navy  with  between  2  and  3  years  of  service.  Similarly  there  were  371 
"Electronics  Technicians"  with  more  than  11  and  not  more  than  12  years  of 
service. 

Using  this  data  and  formula  (6),  we  can  calculate  a  point  estimate  of 
the  continuation  rates  8^  .  For  example  for  Boatswain's  Mates. 


-fff-0.89 

4  n_1  3  382 


on 


Tl0.2  2832 

P2  n  .  .  1547 

-x,i 


1.83  . 


In  this  case,  we  see  that  the  continuation  rate  exceeds  1.  This  is  due  to  a 
large  number  of  inputs  into  BM  rating  with  two  years  of  unrated  service.  The 
continuation  rates  and  survivor  fraction  calculated  for  the  skill  categories 
BM  and  ET  are  presented  in  Tables  3. 

In  the  critical  3kill,  ET,  the  continuation  rates  are  less  than  .1  and 
survivor  fractions  decrease.  This  follows  from  the  Navy's  practice  of  training 
ET's  from  recruits  with  a  zero  LOS,  The  BM  rating  draws  its  input  from  the 
unrated  enlisted  ranks,  and  generally  with  an  LOS  equal  to  2.  The  BM  survivor 
fractions  can  be  modified  by  dividing  them  by  02  =  6.28  .  The  modification 
allows  us  to  interpret  as  the  number  of  accessions  at  time  i  that  even¬ 

tually  enter  the  BM  skill  category,  and  ajx^  as  the  number  in  the  BM  skill 
category  in  period  i  +  j  . 
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TABLE  3 


CONTINUATION  RATES  SURVIVOR  FRACTIONS 


LOS 

BM 

ET 

BM 

ET 

0 

1 

1 

1 

1 

1 

3.43 

0.71 

3.43 

0.71 

2 

1.83 

0.93 

6.28 

0.66 

3 

0.21 

0.85 

1.32 

0.56 

4 

0.89 

0.91 

1.18 

0.51 

5 

0.88 

0.64 

1.04 

0.32 

6 

0.33 

0.85 

0.97 

0.27 

7 

0.95 

0.79 

0.92 

0.22 

8 

0.93 

0.83 

0.85 

0.18 

9 

0.86 

0.90 

0.74 

0.16 

10 

0.96 

0.94 

0.70 

0.15 

11 

0.96 

0.94 

0.68 

0.14 

12 

0.98 

0.99 

0.66 

0.14 

13 

0.95 

0.93 

0.63 

0.13 

14 

0.98 

0.91 

0.62 

0.12 

15 

0.98 

0.99 

0.60 

0.12 

16 

0.96 

0.99 

0.57 

0.12 

17 

0.98 

0.99 

0.57 

0.12 

18 

0.70 

0.79 

0.39 

0.09 

19 

0.52 

0.25 

0.20 

0.02 

20 

0.62 

0.50 

0.12 

0.01 

21 

0.78 

0.99 

0.10 

0.01 

22 

0.65 

0.99 

0.07 

0.01 

23 

0.77 

0.99 

0.05 

0.01 

24 

0.76 

0.99 

0.04 

0.01 

swwssa 
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4.  MINIMIZING  ACCESSION  COSTS 

The  accession  schedule  that  meets  future  manpower  requirements  z exactly 
is  tound  by  solving  the  equations 


Vo 


2i  '  yi 


Vi  +  aox2 


zo  - 


(i) 


VlXl  +  am2X2  +  *  *  ‘  +  Vm 


z  -  y 
m  m 


a  x-  +  a  ,x_  +  ...  +  a.x  +  artx 
m  1  m-1  2  1  m  0  m+1 


m+1 


a  x„  +  . . .  +  a-x  +  a.x +  anx ..  ■  z  .. 
m2  2  m  1  m+1  0  m+2  m+2 


In  this  section  we  shall  relax  the  assumption  that  future  manpower  require¬ 
ments  are  satisfied  exactly.  Instead  we  shall  treat  the  variables  zfc  as 
lower  bounds  on  the  manpower  level  at  time  t  .  Moreover,  we  shall  restrict 
the  accessions  xt  to  be  nonnegative.  Thus  the  equalities  in  (1)  will  be 
replaced  by  inequalities  (>)  .  This  leaves  us  with  an  infinite  system  of 
linear  inequalities  that  will,  in  general,  have  a  large  number  of  possible 
solutions.  To  obtain  a  single  accession  schedule  in  this  case  we  can  specify 
a  performance  criterion  and  then  select  the  accession  schedule  that  optimizes 
that  criterion. 

In  the  analysis  that  follows  we  assume  that  the  performance  criterion  is 
to  minimize  the  present  worth  of  all  future  accession  costs.  This  objective 
is  obtained  by  discounting  future  costs  to  today's  dollars  and  then  summing 
ove-  all  future  periods.  We  obtain  the  extremely  simple  and  useful  result 
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that  the  optimal  accession  schedule  Is  independent  or  the  costs.  In  addition, 
we  will  show  that  there  are  many  conditions  of  practical  interest  when  the 
equality  solution  discussed  in  earlier  sections  is  indeed  an  optimal  solution 
to  the  infinite-horizon  program.  Finally,  we  will  discuss  several  practical 
methods  by  which  the  infinite  program  can  be  truncated  and  approximated  by  a 
finite  program  with  T  rather  than  an  infinite  number  of  periods.  The  method 
of  truncation  that  should  be  used  depends  on  the  conditions  the  user  wishes  to 
assume  for  planning  periods  in  the  distant  future. 

The  Infinite  Horizon  Program 

Thera  are  several  approaches  to  the  solution  of  (1)  when  the  restriction 
of  exactly  meeting  the  mar  _  rer  requirements  is  relaxed.  In  describing  a  cost 
minimization  model  in  which  the  manpower  requirements  are  considered  to  be 
lower  bounds  we  assume  that  manpower  requirements  can  in  fact  be  modelled  by 
such  inequality  restrictions. 

The  problem  is  to  choose  the  nonnegative  vector  (x^,x2,  ...)  that 
satisfies 


Vi 


£  21  "  yl 


Vl  +  “oX2 


±  Z2  -  y2 


xt  >  0 


(2) 


a2x1  +  +  OqX^  j>  Z3  -  y3 


Note  that  it  is  possible  to  consider  xt  as  the  number  of  accessions  above 
some  minimal  level.  Then  z  should  be  interpreted  as  the  requirements  left 
unsatisfied  by  the  minimal  accession  schedule. 

Let  c^  for  k  =  0,1,2 . m  be  the  cost  of  training  and  support  of 
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an  Individual  in  the  lc  year  of  service.  The  discounted  cost  of  an  accession 
is  thus 


c  ‘  '  Vo  +  “lCl6  +  V2S‘'  +  •"  +  V=»{" 


where  6  is  a  discount  factor  less  than  one.  The  total  discounted  cost  of 
an  accession  schedule  (x^,x2,  . ..,)  is  therefore 


w 

(4)  J  cx  5^  ^  *  ex.  +  ficx~  +  6^cx«  +  ... 

j«l  2 

Since  an  accession  in  any  period  j  costs  6^  ^c  ,  it  is  evident  that  the 
value  of  c  will  not  affect  the  optimal  solution.  In  fact,  we  can,  without 
any  loss  in  generality,  use 


00 

l  x  6^  1  -  x  +  fix.  +  S2x  +  ... 
j-1  J  z  j 


as  the  objective  to  be  minimized  subject  to  the  constraints  in  (2). 


The  Dual  Program  and  Optimality  Conditions 


The  dual  of  the  infinite  linear  program  in  (5)  and  (2)  is  to  find  non¬ 
negative  variables  u^,u2,  ...  which 


Maximize  -  y^  +  u2(z2  -  y2)  + 


subject  to  the  inequality  constraints 


ula0  +  u2°l  +  u3a2  +  *  •  *  —  1 


U2a0  +  u3al  +  —  &  ut  Z.  0 


u3aQ  "h  •  •  •  ^ 


(7) 
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A  feasible  solution  of  the  dual  problem  always  exists  and  is  given  by 
t— 1  /  m  '\ 

u  =  y6  ,  where  y  =  l/[  £  a  H  .  This  solution  is  strictly  positive 

\j=0  j  / 

and  satisfies  each  dual  constraint  as  an  equality. 

The  main  results  of  linear  programming  do  not  carry  over  to  infinite 
linear  programs;  the  strong  duality  theorem  can  fail,  [1],  and  in  some  cases 
even  the  weak  duality  then  era  and  the  related  sufficiency  condition  (comple¬ 
mentary  slackness)  can  fail,  [2].  Fortunately,  our  problem  has  a  great  deal 
of  special  structure  and  under  two  reasonable  assumptions  we  shall  establish 
a  useful  lower  bound  on  the  optimal  value  of  the  infinite  problem  and  deter¬ 
mine  that  the  solution  of  (1)  is  optimal  if  it  is  nonnegative,  i.e„  feasible. 


tn  r 

The  second  assumption  is  on  the  m  order  difference  equation  l  a  x  =  z  , 

j=0  3  t 

with  initial  conditions  given  by  (Xq.x^,  . ...  )  •  It:  is  obvious  from 

equation  (5)  in  Section  2  that  the  initial  conditions  are  also  specified  by 
the  legacy  •••»  ym  >  an<*  that  the  unique  solution  of  this  difference 

equation  is  simply  the  solution  of  (1)  in  this  section.  The  difference  equa¬ 
tion  is  called  asymptotically  stable  if  all  the  complex  roots  of  its  charac- 

m  , 

teristic  polynomial  £  a.w  3  ,  have  absolute  value  strictly  less  than  one. 

J-0  3 

In  particular,  the  condition  is  satisfied  if:  1  _>  >_  and  >  “i  >  ai+i 

for  some  i  ,  and  appears  to  be  satisfied  by  other  survivor  fraction  data  that 
we  have  encountered  to  date.  In,  [3],  Chapter  10  it  is  demonstrated  that 
asymptotic  stability  implies  that  for  any  initial  conditions  and  any  bounded 
zt  the  solution  of  the  difference  equation  will  be  bounded.  Let  xt  be  the 
unique,  bounded  solution  of  (1).  From  our  two  assumptions  it  follows  that 
xfc  <=  Max  (O.x^,)  is  a  bounded  and  feasible  solution  for  the  constraints  (2). 

To  establish  optimality  conditions  it  is  useful  to  examine  the  T  period 
approximations  of  the  infinite  problem.  The  T  period  problem  is 
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Minimize  £  6^  *x. 

j-1  3 

t-1 

(8)  Subject  to  l  ax  .  >  z  -  y 

j«0  J  c  J  t  c 

xt  >  0  for  t  -  1,2,  ...,  T 

Notice  that  any  solution  of  the  infinite  problem  will  have  its  first  T  com¬ 
ponents  (x^,x2,  ...,  x^,)  feasible  for  the  T  period  approximation.  Moreover, 
the  solution  ufc  -  pfi1"  *  for  t  ■  1,2,  ...»  T  will  be  feasible  for  the  dual 
of  problem  (8).  Thus  by  the  weak  duality  theorem  for  problem  (8),  we  have 

l  >  l  -  yt) 

t-1  c  t-1 

for  all  T  and  any  solution  of  the  infinite  problem.  This  gives  us  a  lower 

OO 

bound  £  pS*  1(z  -  y  )  on  the  optimal  value  of  the  infinite  problem. 

t-1  ' 

Finally,  since  6  <  1  ,  and  xt  is  bounded  it  is  possible  to  show  that 

00  oo 

£  6*"  Xx  -  £  1(z  -  y  )  .  Thus  the  possibly  infeasible  solution  x 

t-1  C  t-1  1 

achieves  the  lower  bound.  It  follows  that  xt  is  optimal  if  it  is  feat  ible, 

l.e.  if  each  element  of  x^  is  nonnegative. 

In  summary  the  result  of  tnis  section  is  the.t  under  the  hypothesis  of 

asymptotic  stability,  boundedness  of  the  requirements  sequence  zfc  ,  and 

nonnegativity  of  the  survivor  fractions,  the  equality  solution  *t  is  optimal 
if  it  is  nonnegative. 

Optimality  of  Equality  Solution 


Intuition  suggests  the  equality  solution  of  (2)  will  often  be  an  optimal 
solution  of  the  infinite  horizon  planning  problem  posed  in  (5).  Identification 
of  the  conditions  which  must  exist  for  the  equality  solution  to  be  nonnegative 
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and  thus  optimal  gives  considerable  insight  into  the  structure  of  the  problem. 
These  conditions  place  bounds  on  the  magnitude  of  allowable  changes  in  future 
manpower  requirements.  Moreover,  optimality  of  the  equality  solution  gives  us 
an  analytic  expression  for  total  system  cost  as  a  function  of  the  continuation 
rates,  net  requirements,  and  costs. 

To  begin  the  analysis  of  this  section  we  look  for  solution  of  a  reduced 
system  of  linear  equations 


Vi + 

X2 

=  4»1 

Wi + 

61x2  +  X3 

*  *1*2 

B1B2B3X1  + 

ei62X2  +  P1X3  +  X4 

*  <t>1^ 2^ 

which  is  equivalent  to  normalizing  the  equality  system  in  (1)  by  the  constant 
.  In  (9)  <{>.  is  defined  to  be  the  ratio  of  successive  net  require- 


Z1  -  yl 


ments,  i.e.  z.,n  -  y . , ,  /z .  -  y.  .  Alternatively,  one  can  view  (9)  as  the 

J+I  j-t\L  j  j 

equality  system  in  (2)  with  =  1  ,  ^ j  "  aj  *  Z1  "  yl  =  ^  anc* 

-  y  )  =  ^ j  *  e^t^ier  case,  multiplication  of  a  solution 

vector  x  =■  (x^,X2»  •••)  of  (9)  by  z^  -  y^/a^  yields  the  equality  solution 
of  (1).  It  should  be  noted  by  the  reader  that  just  as  the  continuation  rates 
6^,02*  •••  Pj  measure  the  period  to  period  increases  or  decreases  in  the 
fraction  surviving  an  additional  period,  the  growth  rates  •••  ^ 

measure  the  period  to  period  increases  or  decreases  in  the  net  requirements 


that  must  be  met  by  new  accessions. 

With  this  definition  of  terms  it  is  tempting  to  believe  that  so  long  as 
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4> j  >_  Sj  ,  i.e.  the  one  period  increase  (decrease)  in  net  requirements  exceeds 
the  one  period  change  in  the  continuation  rates,  nonnegative  accessions  will 
meet  new  requirements  in  the  next  period.  However,  the  simple  numerical  ex¬ 
ample  having  data 


(10a) 


(♦1.*2.*3.+4)  -  <5  0.2  1  1) 

(81,e2,e3,6,)  *  (2  0.05  1  1) 


yields  a  solution 

¥  • 

(10b)  (x1,x2,x3,x^,x5)  =(13  5.1  10.8  20.49) 


with  negative  components.  This  example  shows  that  the  condition  <J>j  -  8j  0 
is  insufficient  to  guarantee  nonnegative  accessions  in  all  periods  if  there 
are  periods  of  very  large  requirements  followed  by  a  period  of  small  require¬ 
ments. 

It  is  not  difficult  to  show  by  a  direct  substitution  of  unknowns  in  the 
first  three  equations  of  (9)  that  an  equality  solution  x^  satisfies 


(ID 


1 


x3  “  (d>14> 2  "  "  ei^i  “  ei* 


3y  a  reversion  of  the  series  in  (9)  it  can  be  shown  that  in  general 
satisfies  the  j  order  linear,  homogeneous  difference  equation, 


x 


i 


J 

l  X 


j+l  ^  ;)+i-i'rj 


i-l 

(<J>.  -  0,)  n  3 


i'  “  Kk  * 
k=0 


(12) 


j  1  1 
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With  (12)  it  is  now  possible  to  obtain  recursively  in  terms  of 


*1 »X2 ’ 


Xj  .  As  we  have  already  pointed  out  the  solution  of  only 


depends  on  0^'s  and  i^’s  with  i  £  j  .  In  other  words  changing  the 
requirements  on  periods  beyond  j  does  not  affect  the  accessions  in  time 
periods  on  or  before  j  .  Since  continuation  rates  are  nonnegative  it  is 
simple  to  show  that  the  sequence  of  inequalities 


(13)  >  Max  {31,62,  ...  BjJ  j>l 

is  sufficient  to  ensure  that  x^  >_  0  .  The  system  of  inequalities  in  (13) 
is,  of  course,  much  more  restrictive  than  the  one  period  inequalities 

21  Sj  (all  j)  as  it  compares  the  growth  rate  in  one  period  to  all  pre¬ 
vious  continuation  rates.  The  proof  is  straightforward:  x^  =  1  is  non¬ 
negative  and  by  (11)  <j>^  >_  0^  also  implies  x^  >_  0  .  If  we  now  assume  that 

x, ,x_,  ...  x.  are  all  nonnegative  and  that  (13)  holds,  then  x..,  in  (12) 
12  j  j+1 

is  a  sura  of  nonnegative  terms.  Thus,  by  induction  on  n  we  see  that  (13)  is 
sufficient  to  guarantee  that  all  accessions  x^  are  nonnegative. 

Again  we  use  the  data  in  (10a)  to  indicate  why  the  equality  solution  in 
(10b)  fails  to  be  optimal.  Notice  that 


=  5  _>  Max  { 2 }  =  2 

4>2  =  0.2  <  Max  {2,0.05}  =  2 
(l'O 

<f>3  =  1  <  Max  {2,0.05,1}  -  2 
=  1  <  Max  {2,0.05,1,1}  =  2  . 


Necessary  conditions  for  the  in  (12)  to  be 


by  noting  that  we  can  always  rewrite  the  (j  +  1) 


st 


(15) 


J+1 


*1*2*3  "•  *j  -  W3  •••  -  ^  xi 

i=2 


nonnegative  are  obtained 
equation  in  (9)  as 


j+l-i 

n 

k=0 
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Since  the  right  hand  side  is  nonnegative  if  all  x^  _>  0  it  follows  that  the 
inequalities  on  cumulative  products, 

*1±31 

(16)  tj+2  1  6^2 

^1^2^3  -  616283  »  etc* 

must  hold.  While  (fj  0^  impuy  (16)  the  converse  is  not  true  as  it  is  quite 
possible  that  <j>j  <  0^  while  •••  —  3132  3j  *  Thus  8*-mple 

and  local  test  of  whether  the  growth  rate  in  net  requirements  exceeds  the  con¬ 
tinuation  rate  in  each  period  lies  somewhere  between  the  necessary  conditions 
of  (16)  and  the  more  global  sufficiency  conditions  in  (13). 

When  the  equality  solution  is  nonnegative  the  optimal  dual  variables  for 
problem  (7)  are  given  by  u6C  ^  .  However,  to  obtain  useful  total  and  marginal 
cost  information  it  is  necessary  to  reintroduce  (3);,  the  discounted  cost  per 
accession.  The  total  cost  of  an  optimal  solution  is 


(17) 


(f  (y^/5  •/)  I  (..-v**-1 

\j=0  2  1  'j=0  2  /  t=l  c  c 


This  formula  has  a  reasonable  interpretation  in  the  case  where  =  1  ,  and 
the  Oj  are  nonincreasing.  Let  be  the  probability  that  an 

individual's  lifetime,  i.e.  maximum  LOS,  is  equal  to  j  .  With  this  stochastic 
interpretation  of  the  survivor  fractions  we  define  two  random  variables:  T  the 
individual's  lifetime  and  K  the  total  support  cost  of  the  individual.  When 
6  is  equal  to  1,  the  term  in  brackets  in  (17)  is  simply  E[K]/e[T]  .  Thus  we 
can  save  by  keeping  E[T]  fixed  and  reducing  costs.  Notice,  however,  that  if 
we  attempt  to  increase  expected  lifetime  by  changing  the  ,  then  the  cost 
will  change  also.  We  can  get  a  more  accurate  estimate  of  the  impact  of  possible 
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changes  by  rewriting  the  first  terra  of  (17)  with  6=1,  and  expressed 

in  terms  of  continuation  rates 


(m  j  v  m  j 

l  c1  n  8  )/  l  n 

j=0  J  k=0  *7  j=0  k= 


The  derivative  of  (18)  with  respect  to  8  is 

Xj 


m 

l  L.  -  e[K]/e[T]U 
,i=&  '  J  ]  2 


j^rrr 


This  expression  reinforces  some  intuitive  feelings  about  the  system.  If  the 
cost  in  the  periods  following  £  is  greater  than  average,  then  an  increase  in 
8^  will  surely  increase  costs.  On  the  other  hand,  if  all  of  the  downstream 
costs  are  less  than  average,  then  it  will  reduce  costs  when  8^  is  increased. 

Special  Cases 


This  section  examines  several  special  cases  and  derives  tighter  and  more 

easily  verified  conditions  under  which  the  equality  solution  is  optimal. 

First  if  the  8^  are  nondecreasing,  then  4>^  _>  8^  implies  ^  _>  8^  for 

i  <  i  ,  thus  the  conditions  <j>.  >  8  are  necessary  and  sufficient  for  the 

J  ~  j 

equality  solution  to  be  nonnegative. 

In  a  second  case,  if  the  are  nonincreasing,  then  8j  ^  1  for  all  j  . 

Moreo/er,  nonincreasing  imply,  see  equation  5  Section  2,  that  the  legacy 

yt  is  nonincreasing.  If  z t  is  nondecreasing,  it  follows  that  z -  y  is 
nondecreasing  and  thus  that  <J> ^  1  for  all  j  .  Therefore  we  have  optimality 

of  the  equality  solution  under  the  readily  verified  conditions  z^  nondecreasing 
and  nonincreasing. 

If  we  further  specialize  the  first  case  so  that  8^  =  8  <  <*>j  for  all  j  , 
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then  we  can  write 


/  ♦j  -  6  \ 

(19) 

Xj+1  "  -  e ) 

XJ  3  ■ 

2,3,  • . . 

with  ■ 

1  ,  and  x2  -  *1  -  B  .  If 

■  *  j  “  • 

1,2,  • • • , 

thereafter 

we  obtain 

1 

-  e> 

j  <  T 

(20) 

Xj+1  |*T<1  -  8) 

j  >  T 

Calculating  Approximately  Optimal  Accessions  Schedules 

It  is  often  not  possible  to  determine  if  the  equality  solution  of  (1)  is 
nonnegative.  These  situations  typically  Involve  relatively  large  decreases 
in  requirements  in  the  first  few  periods  together  with  a  large  legacy.  In 
these  cases  we  must  find  some  procedure  for  either  solving  or  approximating 
the  solution  of  the  infinite  horizon  problem  (2)  and  (5). 

This  section  will  briefly  describe  three  methods  for  calculating  approxi¬ 
mately  optimal  solution  to  the  Infinite  horizon  optimization  problem.  Each  of 
the  three  methods  is  based  on  a  partition  of  the  original  infinite  problem  into 
a  T  period  finite  problem  followed  by  an  infinite  problem  that  commences  at 
time  T  +  1  .  The  hope  is  that  the  system  will  settle  down  enough  so  that  the 
problem  starting  at  time  T  +  1  will  have  a  nonnegative  equality  solution  re¬ 
gardless  of  the  choice  of  (x^,x2,  ...»  x^  . 

The  first  method  simply  Ignores  the  decisions  and  constraints  for  time 
T  +  1  onwards.  Thus  we  are  solving  problem  (8)  of  this  section.  This  pro¬ 
cedure  is  quite  simple,  however,  it  can  lead  to  optimal  programs  that  save  in 
periods  1  through  T  by  presenting  difficult  initial  conditions  for  the 
second  problem  that  commences  at  time  T  +  1  .  Since  the  problem  that  starts 
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at  time  T  +  1  is  not  explicitly  considered  in  the  objective,  there  is  no 

penalty  to  deter  this  type  of  behavior. 

The  second  method  attempts  to  provide  a  smooth  transition  to  equilibrium 

by  fixing  accessions  at  their  equilibrium  value  for  periods  T  +  1  onward. 

m 

The  assumption  is  that  z  «*  z  for  t  ^  T  ,  and  that  x  =>  z  /  £  a  for 

1  1  t  l  j=0  J 

all  t  >_  T  .  Thus  the  accessions  in  periods  1  -  m,2  -  m,  ....  -1,0  and 
T  +  1,  T  +  2,  ...  are  all  known.  We  must  determine  the  accessions  in  periods 
1  through  T  in  order  to  satisfy  the  lower  bound  requirement  in  the  first 
T  +  m  time  periods.  This  leads  to  a  linear  program  with  T  +  m  inequality 
constraints  and  T  *  nonnegative  variables  x^ . xT  .  The  dual  linear  pro¬ 

gram,  with  T  inequality  constraints  and  T  +  m  nonnegative  variables,  is 
easier  to  solve.  Unfortunately  this  truncation  procedure  has  not  been  effec¬ 
tive  in  numerical  examples  we  have  solved  to  date.  We  frequently  obtain  rela¬ 
tively  low  values  of  x^,_^  ,  xT_2  ,  etc.  and  a  relatively  large  value  of  xT  . 
In  effect,  the  program  satisfies  the  boundary  restriction  by  making  a  last 
period  correction.  This  behavior  is  contrary  to  the  smooth  transition  to 
equilibrium  that  the  model  was  designed  to  produce. 

The  third  method  is  based  on  the  theory  developed  previously  on  optimality 
conditions.  The  lower  bound  on  the  optimal  value  of  the  infinite  horizon  prob¬ 
lem  that  starts  at  time  T  +  1  ,  must  consider  the  additional  legacy  due  to  the 
accessions  in  periods  1  through  T  .  The  bound  is 


CO 


l 

j  =TJr  1 


l 


k=l 


“j-kX 


f  k-1 

If  we  add  the  value  of  the  first  T  periods,  l  6  x,  ,  and  rearrange  terms 

k=l  k 

we  obtain  a  lower  bound  for  the  original  infinite  problem  in  terms  of  the  de¬ 
cision  variables  x^ . x^,  . 
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6  Vi 


nrfk-T 

1 

i- 1 


.t-1 
6  a 


l+T-k)^ 


l 

*T+1 


- 


V 


Notice  that  the  expression  on  the  right  is  a  constant,  independent  of  the 
decisions  in  period  1  through  T  .  Thus  the  linear  program  we  solve  is 

T  /  .  -  nrfk-T  .  \ 

Minimize  j  ( 5  “  -  J  6  "V  wJi 

k-1  '  i-1  K 

t 

(21)  Subject  to  l  at_kxk  >_  z  -  y 

k*l 

x£  ^  0  t  -  1,2,  ...,  T  . 

Notice  that  the  third  method  is  similar  to  the  first,  except  for  the  objective 
which  explicitly  contains  a  penalty  cost  on  the  legacy  created  for  the  infinite 
problem  starting  at  time  T  +  1  . 

In  our  calculations  we  have  used  the  third  method.  There  are  two  main 
reasons  for  this:  first,  the  lower  bound  is  exact  when  the  equality  solution 
ie  optimal  for  the  problem  commencing  in  period  T  +  1  ;  second,  the  calcu¬ 
lations  we  have  made  using  actual  survivor  fractions  have  led  to  realistic 
answers. 

The  dual  linear  program  is  stated  below: 


Maximize  £  u  (z  -  y  ) 
t-1  t  t  C 


T-t 

(22)  Subject  to  £  o  u  +  v 

j-0  J  C  J  C 


,  _  nrft-T  .  . 

fi*"1  -  Jp  l  6s’’"1* 

do  *+t-c 


ufc  0  ,  v  j>  0  for  t  -  1,2 . 


T 
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The  interpretation  of  the  dual  variables  is  straight  forward:  is  the 

marginal  change  in  the  optimal  cost  of  increasing  net  requirements  (zt  -  y  )  , 
and  v  is  the  marginal  change  in  the  optimal  cost  associated  with  increasing 
the  lower  bound  on  accessions  in  period  t  .  The  dual  constraints  can  be  in- 

t-1 


m 


terpreted  as  the  equations  £  =  6 

t  >  T  have  been  assigned  the  values  5  p  and  0  . 


where  ufc  and  vfc  for 


The  example  below  was  solved  for  the  rating  ET,  electronics  technician, 
using  the  survivor  fraction  and  age  distribution  data  contained  in  Tables  2 
and  3  of  Section  3.  We  imposed  a  lower  bound  of  1750  accessions  per  year,  a 
discount  factor  6  =  0.95;  and  z ■  16000  for  t  ^  6  . 


m 

l 

j=0 


1 

2 

3 

4 

5 

zt 

20000 

18000 

16000 

16  000 

16000 

xt 

2112 

1750 

1750 

2098 

2828 

Vt-j 

20000 

18363 

16922 

16000 

16000 

u 

t 

0.5 

0.0 

0.0 

0.19 

0.18 

vt 

0.0 

0.35 

0.2 

0.0 

0.0 

Since  z  is  constant  for  t  >  6  and  the  survivor  fractions  a.  are 
t  —  j 

nonincreasing  for  the  ET  rating  we  know  that  the  equality  solution  is  optimal 
for  the  problem  beginning  in  period  6  regardless  of  the  accessions  in  periods 
1  through  5.  Thus  the  solution  above  is  optimal  for  the  infinite  horizon  problem. 


I 

It 

I 

fife 
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